Search for: All records

In weathered bedrock aquifers, groundwater is stored in pores and fractures that open as rocks are exhumed and minerals interact with meteoric fluids. Little is known about this storage because geochemical and geophysical observations are limited to pits, boreholes, or outcrops or to inferences based on indirect measurements between these sites. We trained a rock physics model to borehole observations in a well-constrained ridge and valley landscape and then interpreted spatial variations in seismic refraction velocities. We discovered that P-wave velocities track where a porosity-generating reaction initiates in shale in three boreholes across the landscape. Specifically, velocities of 2.7 ± 0.2 km/s correspond with growth of porosity from dissolution of chlorite, the most reactive of the abundant minerals in the shale. In addition, sonic velocities are consistent with the presence of gas bubbles beneath the water table under valley and ridge. We attribute this gas largely to CO₂produced by 1) microbial respiration in soils as meteoric waters recharge into the subsurface and 2) the coupled carbonate dissolution and pyrite oxidation at depth in the rock. Bubbles may nucleate below the water table because waters depressurize as they flow from ridge to valley and because pores have dilated as the deep rock has been exhumed by erosion. Many of these observations are likely to also describe the weathering and flow path patterns in other headwater landscapes. Such combined geophysical and geochemical observations will help constrain models predicting flow, storage, and reaction of groundwater in bedrock systems.

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east‐west trending shale watershed in the Valley and Ridge province in the USA Using wave equation travel‐time tomography from a seismic array of >4,000 geophones, we obtained a 3D P‐wave velocity (Vp) model that resolves structures ∼20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2,700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2,700 m/s contour is greater under the N‐ as compared to S‐facing hillslopes and under sub‐planar as compared to concave‐up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N‐ versus S‐facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P‐wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. The Vp result is best explained as the long‐term integrated effect of groundwater flow‐induced geochemical weathering of shale in response to climate‐driven patterns of micro‐ and macro‐topography.

We examine upper mantle anisotropy across the Antarctic continent using 102 new shear wave splitting measurements obtained from teleseismic SKS, SKKS, and PKS phases combined with 107 previously published results. For the new measurements, an eigenvalue technique is used to estimate the fast polarization direction and delay time for each phase arrival, and high‐quality measurements are stacked to determine the best‐fit splitting parameters at each seismic station. The ensemble of splitting measurements shows largely NE‐SW‐oriented fast polarization directions across Antarctica, with a broadly clockwise rotation in polarization directions evident moving from west to east across the continent. Although the first‐order pattern of NE‐SW‐oriented polarization directions is suggestive of a single plate‐wide source of anisotropy, we argue the observed pattern of anisotropy more likely arises from regionally variable contributions of both lithospheric and sub‐lithospheric mantle sources. Anisotropy observed in the interior of East Antarctica, a region underlain by thick lithosphere, can be attributed to relict fabrics associated with Precambrian tectonism. In contrast, anisotropy observed in coastal East Antarctica, the Transantarctic Mountains (TAM), and across much of West Antarctica likely reflects both lithospheric and sub‐lithospheric mantle fabrics. While sub‐lithospheric mantle fabrics are best associated with either plate motion‐induced asthenospheric flow or small‐scale convection, lithospheric mantle fabrics in coastal East Antarctica, the TAM, and West Antarctica generally reflect Jurassic—Cenozoic tectonic activity.

We used seismic refraction to image the P‐wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill‐valley‐hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. This contradiction could be explained by the greater importance of infiltration‐driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.